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We study a 1-dimensional AKLT spin chain, consisting of spins S in the bulk and 
S/2 at both ends. The unique ground state of this AKLT model is described by the 
Valence-Bond-Solid (VBS) state. We investigate the density matrix of a contiguous 
block of bulk spins in this ground state. It is shown that the density matrix is a 
projector onto a subspace of dimension (S + l) 2 . This subspace is described by non- 
zero eigenvalues and corresponding eigenvectors of the density matrix. We prove 
that for large block the von Neumann entropy coincides with Renyi entropy and is 
equal to In (S + l) 2 . 
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1. Introduction 

There is considerable current interest in studying various interacting quantum sys- 
tems from the quantum information perspective. Quantum entanglement is a fun- 
damental measure of how much quantum effects we can observe and use to control 
one quantum system by another, and it is the primary resource in quantum com- 
putation and quantum information processing (Bennett & DiVincenzo 2000, Lloyd 
1993). Entanglement properties also play an important role in condensed matter 
physics, such as phase transitions (Osterloh, et al. 2002; Osborne & Nielsen 2002) 
and macroscopic properties of solids (Ghosh, et al. 2003; Vedral 2004). Extensive 
research has been undertaken to understand quantum entanglement for spin chains, 
correlated electrons, interacting bosons as well as other models, see Amico, et al. 
(2007), Audenaert, et al. (2002), Fan & Korepin (2008), Katsura, et al. (20076), 
Fan, et al. (2007), Arnesen, et al. (2001), Korepin (2004), Verstraete, et al. (2004a, 
b), Campos Venuti, et al. (2006), Jin & Korepin (2004), Vedral (2004), Latorre, et 
al. (2004a, b, 2005), Orus (2005), Orus & Latorre (2004), Pachos & Plenio (2004), 
Plenio et al. (2004), Fan & Lloyd (2005), Chen, et al. (2004), Zanardi & Rasetti 
(1999), Popkov & Salerno (2004), Keating & Mezzadri (2004), Gu, et al. (2003, 
2004), Wang, et al. (2004), Wang & Kais (2004), Holzhey, et al. (1994), Calabrese 
& Cardy (2004), Levin & Wen (2006), Kitaev & Preskill (2006), Ryu & Hatsugai 
(2006), Hirano & Hatsugai (2007) for reviews and references. Characteristic func- 
tions of quantum entanglement, such as von Neumann entropy and Renyi entropy, 
are obtained and discussed through studying reduced density matrices of subsys- 
tems (Fan, et al. 2004; Katsura, et al. 2007a; Its, et al. 2005; Franchini, et al. 2007, 
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2008; Vidal, et al. 2003). An area law for the von Neumann entropy in harmonic 
lattice systems has been extensively studied (Cramer, et al. 2006, 2007; Plcnio, et 
al. 2005). 

In this paper we study a spin chain model introduced by Affleck, Kennedy, Lieb 
and Tasaki (AKLT) (Affleck, et al. 1987, 1988). We consider a 1-dimensional AKLT 
model with spin-S* in the bulk and spin-S*/2 at both ends. The ground state of this 
model is a unique pure state (Arovas, et al. 1988). It is known as the Valence-Bond- 
Solid (VBS) state, which plays a significant role in condensed matter physics. The 
VBS state is closely related to Laughlin ansatz (Laughlin 1983, Iblisdir et al. 2007) 
and fractional quantum Hall effect (Arovas, et al. 1988). It enables us to understand 
ground state properties of anti-ferromagnetic integer-spin chains where the finite 
energy gap known as the Haldane gap exists (Haldane 1983). Universal quantum 
computation based on VBS states (Verstraete & Cirac 2004) and an implementation 
of the AKLT hamiltonian in optical lattices (Garcia-Ripoll et al. 2004) have also 
been proposed. 

The density matrix of a contiguous block of bulk spins (we call it the density 
matrix later for short) has been studied extensively in Kirillov & Korepin (1990), 
Verstraete, et al. (2004a), Fan, et al. (2004), Katsura, et al. (2007a), Freitag & 
Mullcr-Hartmann (1991). It contains information of all correlation functions (Jin 
& Korepin 2004, Katsura, et al. 2007a, Arovas, et al. 1988). Moreover, it has been 
shown in Fan, et al. (2004), Katsura, et al. (2007a) that the density matrix is 
independent of the size of the chain and the location of the block relative to the 
ends. Therefore we can take the length of the block equal to the length of the whole 
chain, (i.e. we can add two ending spins S/2 directly to the block.) Then by using 
the Schmidt decomposition (Nielsen & Chuang 2000), we can show that the density 
matrix of the block is equivalent to the density matrix of the two ending spins. By 
equivalent we mean that all non-zero eigenvalues are the same. Using this method, 
eigenvalues of the density matrix as well as entanglement entropies were obtained 
(Fan, et al. 2004; Katsura, et al. 2007a; Freitag & Muller-Hartmann 1991) without 
knowing the eigenvectors explicitly. 

However, eigenvectors of the density matrix have their own importance. They 
can be used to study the structure and symmetries of the density matrix explicitly 
both for finite block and in large block limit. The construction of eigenvectors also 
provides us with a possible method to diagonalize the density matrix directly. As 
to be shown in following sections (see §2c, §2d, §4c), the eigenvectors also have 
their own physical meaning as degenerate zero-energy ground states. In the context 
of the Haldane gap, these degenerate states are known as edge states and have been 
observed in the S = 1 spin chain compound (Hagiwara et al. 1990). Furthermore, 
eigenvectors become indispensable in quantum computing algorithms, particularly 
in discussing quantum measurements. 

In this paper, we consider AKLT models with two different boundary conditions. 
Let's first take spin 5 = 1 for example. The system consists of a linear chain of N 
spin-l's in the bulk, and two spin-l/2's on the boundaries. We shall denote by Sj 
the vector spin operator at site j (j = 0, 1, . . . , N + 1). The Hamiltonian is 




(1.1) 
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The boundary terms ir describe interactions of a spin- 1/2 and a spin-1. Each term 
is a projector onto a state with spin 3/2: 

2 2 

7T0,1 = - (1 + Sq ■ Si) , 7T.Ar,jV+l = g (1 + 5jV ' S>jV+l) • (1.2) 

The Hamiltonian (1.1) has a unique ground state (VBS state), thus we shall call it 
the unique Hamiltonian. Alternatively, if we consider spin-l's at every site including 
the boundaries, then the Hamiltonian takes the form 




The ground states of this Hamiltonian are 4-fold degenerate. We shall call (1.3) the 
degenerate Hamiltonian. 

For generic spin-5, the unique Hamiltonian is 

JV-l 2S 
H-uniq ^ ^ ^ ^ C j 

P 3,j+1 + ^0,1 + kn,n+i, (1-4) 

j=l J=S+1 

where the projector P/j+i projects the bond spin «7j,j+i = Sj + Sj+i onto the 
subspace with total spin J (J = 5 + 1, ...,25). The boundary terms describe 
interactions between a spin-5/2 and a spin-5: 

3S/2 35/2 

TT0,1= Yl D J P 0,1> ^N,N+1= D J P N,N+1- ( L5 ) 

J=5/2+l J=5/2+l 

Both coefficients Cj and Dj can take arbitrary positive values. Correspondingly, 
the degenerate Hamiltonian with spin-5 at every site takes the form 

N-l 2S 
j=l J=S+1 

The degeneracy of the ground states is (5 + l) 2 . This will be important in descrip- 
tion of eigenvectors of the density matrix (see %2d and §3). 

Consider the AKLT spin chain system with the unique Hamiltonian (1.4) in the 
VBS ground state. The density matrix p of the whole chain is a projector onto the 
unique VBS ground state (see (2.2)). If we pick up a block of L contiguous bulk 
spins as a subsystem and trace out all degrees of freedom outside the block, then we 
obtain the density matrix p L of the subsystem (see (2.6)). Because of entanglement 
with spins outside the block, p L will no longer be a pure state density matrix 
as p is in general. We shall prove that the density matrix p L is a projector onto 
a (5 + l) 2 -dimensional subspace of the complete Hilbert space associated with 
the block (see §2<i and §3). The degenerate Hamiltonian (1.6) becomes essential 
in description of this subspace. When the degenerate Hamiltonian has its size N 
equal to that of the block L, it is referred to as the block Hamiltonian and denoted 
by Ht which is defined by (2.13). It turns out that the block Hamiltonian H^ (i.e. 
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the degenerate Hamiltonian Hd eg in (1.6) with N = L) defines the density matrix 
p L completely in the large block limit L — > oo. The zero-energy ground states of 
the block Hamiltonian Hi, span the subspace that the density matrix p L projects 
onto. So that p L can be represented as the zero-temperature limit of the canonical 
ensemble density matrix defined by H b : 

e -0H b 

p T = lim ——, — wttt, L — > oo, (1.7) 

where 

L-l 25 

H b = H deg (withN = L) = Y / E Cl ■ ■ ( L8 ) 

j=l J=S+1 

In the zero-temperature limit, contributions from excited states of H b all vanish 
and the right hand side of (1.7) turns into a projector onto the ground states of the 
block Hamiltonian. 

As main subjects of the paper, we will construct eigenvectors and derive expres- 
sions for corresponding eigenvalues of the density matrix. We will show that the 
density matrix is a projector. The paper is divided into four parts: 

1. We calculate the density matrix, prove a theorem on eigenvectors and express 
eigenvalues in two different forms using the Schwinger representation (§2). 

2. We investigate the structure of the density matrix in the large block limit. As 
characteristic functions of quantum entanglement, the von Neumann entropy 
and the Renyi entropy are obtained in the limit (§3). 

3. We study the density matrix using a different representation (a pure algebraic 
method) for spin S = 1 (§4). 

4. An alternative proof of the theorem on eigenvectors is given as we take a 
different approach (§5). 



2. Density Matrix for Generic Spin-S 1 

(a) Ground State of the Unique Hamiltonian 

We start with the ground state of the unique Hamiltonian (1.4). It is given in 
the Schwinger representation by the VBS state (Arovas, et al. 1988) 

N g 

|VBS> = J] - bWi+i) l vac >' ( 2J ) 

3=0 

where a) , are bosonic creation operators and |vac) is destroyed by any of the 
annihilation operators a, b. These operators satisfy [aj,aj] = = Sij with all 

other commutators vanishing. The spin operators are represented as = a/jbj, 

Sj = bjdj, Sj = (Ojdj — bjbj)/2. To reproduce the dimension of the spin-S 1 Hilbert 
space at each site, an additional constraint on the total boson occupation number 
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is required, namely (a]aj + bjbj)/2 = 5. More details and properties of the VBS 
state in the Schwinger representation can be found in Kirillov & Korepin (1990), 
Arovas, et al. (1988), Auerbach (1998). The pure state density matrix of the VBS 
ground state (2.1) is 

p= ]VBSKVBS[_ (2 . 2) 
P (VBS|VBS) 1 ; 

For normalization (VBS|VBS) of the VBS state, see Appendix A. 



(b) Density Matrix of a Block of Bulk Spins 

We take a block of L contiguous bulk spins as a subsystem. In order to calculate 
the density matrix of the block, it is convenient to introduce a spin coherent state 
representation. We introduce spinor coordinates 

{9-4. 9 -A 
(u,v) = f cos -e 1 2, sin -e" 1 2 J , < 9 < w, < <p < 2n. (2.3) 

Then for a point ft = (sin 9 cos <f>, sin 9 sin <j>, cos 9) on the unit sphere, the spin-5 
coherent state is defined as 

m s 7m <2 - 4) 

Here we have fixed the overall phase (a U(l) gauge degree of freedom) since it has 
no physical content. Note that (2.4) is covariant under SU(2) transforms (see §3). 
The set of coherent states is complete (but not orthogonal) such that (Freitag & 
Muller-Hartmann 1991; Arecchi, et al. 1972) 



25 ' ! " 



^-JdCl\Vl)(Cl\^ ]T \S,m)(S,m\=I 2S+1 , 



(2.5) 



where \S, m) denote the eigenstate of S 2 and S z , and hs+i is the identity of the 
(25 + l)-dimensional Hilbcrt space for spin-5. The completeness relation (2.5) can 
be used in taking trace of an arbitrary operator. 

Now we calculate the density matrix of a block of L contiguous bulk spins in 
the VBS state (2.1). By definition, this is achieved by taking the pure state density 
matrix (2.2) and tracing out all spin degrees of freedom outside the block: 

Pl = Tr ,i,..., k -i,k+L,...,N,N+i [p] , l<k, k + L-l<N. (2.6) 

Here the block of length L starts from site k and ends at site k + L — 1. p L 
is no longer a pure state density matrix because of entanglement of the block 
with the environment (sites outside the block of the spin chain). It was shown in 
Section 2 of Jin & Korepin (2004) that entries of the density matrix are multi-point 
correlation functions in the ground state. The original proof was for spin 5 = 1/2. 
This statement is generalized to generic spin-5 in Appendix D. 
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Using the coherent state representation (2.4) and completeness relation (2.5), 
p L can be written as (Katsura, et al. 2007a) 



Pl 



(2.7) 



/ 


k-1 N+1 

n n d ^ 

J=0 j=k+L 


fc-2 N 

n n 

j=0 j=k+L 


i(l-^..O j+1 ) 


5 

B t |VBS L )(VBS L |S 




"(25+1) 
A-k 


'I J 


f 


N+1 


TV 

n 


i(l-%.fj j+1 ) 


s 



Here the boundary operator B and block VBS state |VBSl) are defined as 

S S 

B = (u k -ib k — Vk-idk) {cik+L-iVk+L — bk+L-iUk+h) , 



k+L-2 

ivbs L )^ n 

j=k 



°5 & 5+i- & 5°5+i) i vac )' 



(2.8) 
(2.9) 



respectively. Note that both B and |VBSl) are SU(2) covariant (see §3). The 
expression (2.7) can be simplified. We can perform the integrals over Clj (j = 
0, 1, . . . , k — 2, k + L + 1, . . . , N, N + 1) in the numerator and all integrals in the 
denominator (see Appendix A). After integrating over these variables, the density 
matrix p L turns out to be independent of both the starting site k and the total 
length L of the block. This property has been proved in Fan, et al. (2004) for spin 
5 = 1 (using a different representation, namely the maximally entangled states, see 
§4) and generalized in Katsura, et al. (2007a) for generic spin-5. Therefore, we can 
choose k = 1 and the density matrix takes the form 



Pl 



5 + 1 



1 L 



(25+1)! 



/d^odA L +iS t |VBS L )(VBS L |B (2.10) 

71") J 



(5_+l) 

(4 



with 



B f = (u*b\ - v*a\^ (a\v* L+1 - b[u* L+1 ^ , 
|VBS i )=n(aj6j +1 -6jaJ +1 ) S |vac). 



(2.11) 
(2.12) 



i=i 



The last two integral of (2.10) can be performed, but we keep its present form for 
later use. 



(c) Ground States of the Block Hamiltonian 

In order to describe the eigenvectors of the density matrix (2.10), we first study 
the zero-energy ground states of the degenerate Hamiltonian defined in (1.6). We 
choose the length of the spin chain equal to that of the block, i.e. N = L, then the 
degenerate Hamiltonian is called the block Hamiltonian and reads 

L-l 25 

H b = H deg (withN = L) = J2 J2 C J P ?J+i- ( 2 - 13 ) 

j=l J=S+1 
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Now we define a set of S + 1 operators covariant under SU (2) 

A\ = (ua\ + vb\Y (a\b\ - b\a{^ (ua\ + vb^ , < J < S. (2.14) 

These operators act on the direct product of Hilbert spaces of spins at site 1 and 
site L. Then the set of ground states of (2.13) can be chosen as 

|G; J, ft) = A^VBSl), J = 0,...,S. (2.15) 

Any state |G; J, ft) of this set for fixed J and ft is a zero-energy ground state of 

(2.13) . To prove this we need only to verify: (i) the total power of a\ and b\ is 2S, 
so that we have spin-S 1 at the first site; (ii) — S < Jf 2 = Sf + < S by a binomial 
expansion, so that the maximum value of the bond spin Ji i2 is S (from SU(2) 
invariance, see Arovas, et al. 1988). These properties are true for any other site j 
and bond + 1), respectively. Therefore, the state |G; J, ft) defined in (2.15) has 
spin-S 1 at each site and no projection onto the Jjj+i > S subspace for any bond. 

The set of states {|G; J, ft)} depend on a discrete parameter J as well as a 
continuous unit vector ft. States with the same J value are not orthogonal. The 
rank of a set of states with the same J value is 2J + 1, which can be obtained from 
the completeness relation (B 7) (see Appendix B and Hamermesh 1989). Thus the 
total number of linearly independent states of the set {|G; J, ft)} is ^2j =0 (2J + 1) = 
(S + l) 2 , which is exactly the degeneracy of the ground states of (2.13). So that 
{|G; J, ft)} forms a complete set of zero-energy ground states. 

We also introduce an orthogonal basis in description of the degenerate zero- 
energy ground states. It is shown in Appendix B and Hamermesh (1989) that A^j 

(2.14) can be expanded in terms of spin creation operators (M = — J, . . . , J) 
defined in (B5). Operator ^ JM acts on the direct product of two Hilbert spaces of 
spins at site 1 and site L (B 6) and can be expressed in terms of bosonic creation 
operators in the Schwinger representation (Bl). If we define a set of degenerate 
VBS states {|VBS L (J, M))} such that 

|VBSi(J,M)) = *^|VBS L ), J = 0,...,S, M = -J,...,J, (2.16) 

then these (S+l) 2 states (2.16) are not only linearly independent but also mutually 
orthogonal (Appendix C). Furthermore, any ground state |G; J, ft) can be written 
as a linear superposition over these degenerate VBS states, and vice versa (see (B 8) 
of Appendix B). The set {|VBS L (J, M))} differs from {|G; J, ft)} by a change of 
basis, so that it also forms a complete set of zero-energy ground states. 

(d) Eigenvectors of the Density Matrix 

Eigenvalues of the density matrix (2.10) are derived for spin-1 in Fan, et al. 
(2004) and for spin-5 in Katsura, et al. (2007a). Because the density matrix is 
independent of both the total length of the spin chain and the starting site of the 
block, we can add boundary spins directly to the ends of the block. It was shown 
in Fan, et al. (2004), Katsura, et al. (2007a) by a Schmidt decomposition (Nielsen 
& Chuang 2000) that non-zero eigenvalues of the density matrix (2.10) are equal to 
those of the density matrix of the two boundary spins. All other eigenvalues of the 
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density matrix (2.10) are zero. This fact reveals the structure of the density matrix 
as a projector onto a subspace of dimension (S + l) 2 . 

Now we propose a theorem on the eigenvectors of the density matrix given by 
(2.10). The explicit construction of eigenvectors allows us to diagonalize the density 
matrix directly. The set of eigenvectors also spans the subspace that the density 
matrix projects onto. 

Theorem 2.1. Eigenvectors of the density matrix p L (2.10) with non-zero eigen- 
values are given by the set {\G; J, ft)} (2.15), or, equivalently, by the set 
{\ VBSl(J, M)}} (2.16). i.e. they are zero-energy ground states of the block Hamil- 
tonian H b (2.13). 

We prove the theorem by showing that the density matrix p L (2.10) can be 
written as a projector in diagonal form onto the orthogonal degenerate VBS states 
{|VBSl( J, M))} introduced in (2.16). An alternative proof taking a different ap- 
proach is given in §5. 

First, it is realized from the definition of spinor coordinates (2.3) that if we 
change variables (u, v) to (iv* , — hi*), then the unit vector Q is inverted about the 
origin to — Cl. So that we have (Katsura, et al. 2007a) 

(u*&t - uV) s |vac) = i 5 \/5! | - Cl), (2.17) 

where | — Cl) means a spin-S'/2 coherent state for a point opposite to Cl on the unit 
sphere. Therefore, taking expressions of the boundary operator (2.11) and the 
block VBS state |VBS L ) (2.12), we have 

B^VBSl) = (2.18) 

L-l s 

Sl II { a ] b l+i - b l a ]+i) I ~ ^o) 1 ® l vac ) 2 ® ■ ■ ■ ® |vac) L _i <g> | - 
Consequently the density matrix p L (2.10) can be re- written as 



Pl 



5 + 1 



(25+1)! 



S\S\ 



3 = 1 



L-l 

■■4+1 <S> |vac) 2 (vac| ® • • • <g> |vac)i_i(vac| ® 4+i Y\ ( a jbj+i - bjdj + i) S , 

where 4+i and 4+i are (S + l)-dimensional identities associated with site 1 and 
site L, respectively. In obtaining (2.19), we have changed integral variables from Cl 
, Ql+i to —Clo, — &l+i and performed these two integrals using the completeness 
relation (2.5). Next we notice that (see Appendix B) 

s J 

4+>4+i EE \J,M) h L(J,M\ (2.20) 

J=0M=-,7 
s J 

= E E *JMl vac )i( vac l ® | vac) L (vac |* jm- 
j=om=-j 
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As a result, combining (2.19) and (2.20), recalling definitions of |VBSl) (2.12) and 
|VBSl(J, M)) (2.16), the density matrix p L takes the following final form 



Pl 




L s\s\ s 



?E E * t J M|VBS i )(VBS L |* JM (2.21) 



s\s\ 



J=0 M=-J 
S J 



?E E |VBS L (J,M))(VBS L (J,M)|. 



J=0 M=-J 

The set of degenerate VBS states {|VBS L (J, M)» with J = 0, . . . , 5 and M = 
— J, . . . , J forms an orthogonal basis (see Appendix C). These (5 + l) 2 states also 
forms a complete set of zero-energy ground states of the block Hamiltonian (2.13) 
(see §2c and §Appendix B). So that in expression (2.21) we have put the density 
matrix as a projector in diagonal form over an orthogonal basis. Each degenerate 
VBS state |VBSl(J, M)) is an eigenvector of the density matrix, so as any of the 
state |G; J, Ct) (because of the degeneracy of corresponding eigenvalues of the density 
matrix, see §2 e and §2/). Thus we have proved theorem 2.1. 



(e) Eigenvalues of the Density Matrix (Recurrence Formula) 

Having constructed eigenvectors, we need to specify the corresponding eigenval- 
ues. An explicit expression of eigenvalues is obtained in §2 / . In this subsection we 
express eigenvalues through a conjectured recurrence formula as in Katsura, et al. 
(2007a) and Freitag & Muller-Hartmann (1991). Let's apply the density matrix p L 
(2.10) to the state |G; J,0) (2.15) and get 



p L \G;J,h) 
5+1 



(2.22) 



(25+1)! 



5 + 1 

(4tt) 2 



df2 d^L+iB t |VBS L )(VBS L |BA t J |VBS L ). 



Using the coherent state representation (2.4) and completeness relation (2.5), the 
factor (VBS l |.B,4^|VBSl) in (2.22) can be re-written as 



(VBS L |fl4|VBS L > 



(25+1)! 



(2.23) 



47T 



L-l r 




i S 



-(i-JVfi j+ i) 



(u Vi - VqUiY 



■ (uu\ + vvl) J {u\v* L - vlu* L ) S J (uu* L + vv* L ) J (u L v L+ i - V L U L+ lf . 
r -iS 

The factor |(1 — flj ■ fij+i) under the integral of (2.23) can be expanded in 

terms of Legendre polynomials and further in terms of spherical harmonics as (Kat- 
sura, et al. 2007a; Freitag & Muller-Hartmann 1991) 



■(i-n j -n j+1 ) 



5 + 

47T 



1 5 

— £(2/ + l)A(/,S)fl(JVfi j+1 ) 



1=0 

s 



5 + 



t ]Ta(z,5) J2 Y lm (nj)Y* n (n J+1 ) (2.24) 



z=o 
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with coefficients X(l, S) given by 

_ (-1)'5!(5 + 1)! 

x{l > s) = (s-iy.(s + i + iy.- (2 ' 25) 

Using the expansion (2.24) and orthogonality of spherical harmonics, the integrals 
over Clj with j = 2, . . . , L — 1 in (2.23) can be performed. The result is 



S+l 



L S 



(VBS L |i?4|VBS L ) = ^Ti ^±121 ^(2/ + l)A i - 1 G,5) (2.26) 



/ 



dClidCl L Pi(Cli ■ Cl L ) (u Vi — v Ui) S (uul + vvl) J 



{u\v* L - v\u* L ) S J {uul + K) J {ulVl+i - v L u L+1 ) S . 

We plug the expression (2.26) into (2.22). Using transformation properties under 
SU(2) and binomial expansion (see §3), the integral over f^o yields that 

J d&o (u* b\ - Vo»t) ( u oVi - woUi) S = ^ ^ (uia\ + vib\^j (2.27) 

Similarly we can perform the integral over As a result, the following expres- 

sion is obtained from (2.22): 

1 5 

p L \G; J, 0) = — ^ £(2/ + l)A i_1 (/, S)Kj(Cl) |VBS L > . (2.28) 

^ > 1=0 

The operator Kj(Ct) involved in (2.28) is defined as 

Kj(Q) = J dfiidftz, (u ia \ + Vl b\y (uul + ««i) J - "I«1) S_J 

• (ww* + vv* L ) J (u L a[ + v L bty S P^ ■ Cl L ). (2.29) 

It is expressed as an integral depending on the order I of the Legendre polynomial 
Pi(&i-(Il). kJ((i) can be calculated from the lowest few orders (see §3 for example). 
It becomes increasingly difficult to perform the integral as order / increases. Based 
on the eigenvalues of the density matrix obtained in Fan, et al. (2004), Katsura, et 
al. (2007a), we make a conjecture on the explicit form of the operator Kj(Cl) for 
generic order /: 
Conjecture 1 

K\{h) = (5^1) 2 * (\j(J + 1) - \S(\S + 1)) 4. (2.30) 
Here the polynomial /; (x) satisfy the recurrence relation 

T / \ 2Z + 1 (Ax ,\ T , . I fS-l + l\ 2 T , N ,„„.,, 
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with J = 1 and 7i = ( |* 1)2 • 

Note that it is important that Jf/(0) ex defined in (2.14) and Ii(x) has the 
same order as the Legendre polynomial P/(x). The recurrence relation (2.31) was 
proposed in Freitag & Muller-Hartmann (1991) and used in Katsura, et al. (2007a) 
to obtain the eigenvalues of the density matrix. (The original definition of Ii{x) 
differed from our definition in (2.31) by a factor of ^$r-) Conjecture 1 (2.30) is an 
alternative form of theorem 2.1, which also gives eigenvalues through the recurrence 
relation (2.31). Indeed, expressions (2.28), altogether with (2.30) and (2.31) yields 
that 

p L \G;J,h) (2.32) 

= jsTTf i {21 + 1)A ^ 1(/ ' s)Tl Q J(J + l s{ l s + 1} ) |G; J ' h) - 

Non-zero eigenvalues ( J — 0, 1, . . . , S) are seen from (2.32) as 

A (J) = 7^)2 X> + i)^" 1 ^' S ^ (l J ( J + !) - l S (\ S + !)) • ( 2 - 33 ) 

Since all other eigenvalues of the density matrix are vanishing, then we conclude 
again that the density matrix p L (2.10) is a projector onto a subspacc of dimension 
(S + l) 2 . This subspace is spanned by the set of vectors {|G; J, tl)} (2.15). (The 
rank of the set is equal to (S + l) 2 .) Furthermore, we observe from (2.33) that 
non-zero eigenvalues A(J) depend only on J, not on £l. Therefore, {|G; J, O)} with 
fixed J value spans a degenerate subspace with the same eigenvalue. 



(/) Eigenvalues of the Density Matrix (Normalization of Degenerate VBS States) 

Based on the diagonalizcd form (2.21), eigenvalues of the density matrix p L can 
be derived from the normalization of degenerate VBS states. We obtain an explicit 
expression for eigenvalues in terms of Wigner 3j-symbols in this subsection. 

First, the following property is important: Normalization of the degenerate VBS 
state |VBSl(J, M)) depends only on J and is independent of M. With the intro- 
duction of total spin operators of the block S^ t , S^ ot and Sf ot (see Appendix C), 
we prove the statement as follows: 

(VBS L (J, M ± 1) | VBS L (J, M ± 1)) 

- (J T M)(J ± M + 1) < VBS ^ J ' M)|fi£A|VBSt(J, M)> 

= (J T M)(J±M+1) (VBSl(J ' M)l(St2ot " ( ^ ot)2 T ^tot)|VBS L (J, M)) 
= (VBS L (J,M)|VBS L (J,M)>. (2.34) 

Here we have used the fact that |VBSz,(J, M)) is the eigenstate of S^ ot and S^ ot 
with eigenvalues J(J + 1) and M, respectively (see Appendix C). 

It is also realized that normalization of |VBSl(J, M)) can be calculated from 
integrating the inner product of |G; J, 0,} with itself over the unit vector f2 such 
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that 



- / dCl{G; J,Q\G; J, CI) 
or J 



4ir 

(S + J + !)!(£ - J)\J\J\ 
(2J+1)! 



(2.35) 



(VBS L (J,M)|VBS L (J,M)). 



In obtaining this relation (2.35) we have used expansion (B8) and orthogonality 
(B2) in Appendix B. 

Let's consider the integral involved in (2.35). Using coherent state representation 
(2.4) and completeness relation (2.5) as before, we obtain 



47T 

1 

47T 



Q\G;J,Sl) 



(2.36) 



(25+1) 



47T 



(1-SV^l) 



dil 



S-J r 



L 



-(l + Oi - Q) 



L-l 

n 



■(i-Ci r Ci j+1 ) 



-a + Ci-Ci L ) 



Now we expand 



|(1 — Qi ■ Clj) in terms of spherical harmonics as in (2.24), then 
integrate over Cl and from Cl 2 to the right hand side of (2.36) is equal to 

S S-J J h l L I 



4tt((2S'+1)!) l 



1)2 51 



^l— ^z,— /— mi- — /i mL — — lji m— — l 
,L-1/; cAW; C T\\2/ 



(S + 1) L -!(S- J+1)(J + 1) 2 

• y dfii J dCl L \ L - 1 (l 1 ,S)\(l L ,S - J)\\l, J) 
■ Y lu ^{Sli)Y lLtmL <Slx)Y l>m <Sl^ (2.37) 
Here we apply the following useful formula: 

J dClY lumi (Q)Y lLtmL (n)Y ltm (Cl) 

'97-, 4- 1V97, 4- 1V97-I- 1 W 7, 7. / \ / 7, K 7 \ 

,(2.38) 



(2ii + l)(2/i + l)(2/ + l) f Z! f L i 





47T 



mi m l m 



where 



k h I 
mi rah m 
out the integrals in (2.37) and obtain 



is the Wigner 3j'-symbol. Using formula (2.38), we carry 
and obt 

((2S+1)\) L 



s S-J J h 



(s + i^-hs- j + i)(j + i) 



1)2 51 



/l=0 £z,=0 ^— mi- — Zi m L — — l L m— — l 

(2k + 1)(2Z L + 1)(2/ + l)A i_1 (Zi, S)\(Il, S - J)A 2 (Z, J) 

2 / , , , x 2 







Ji k I 
mi mi m 



(2.39) 
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The symbols obey the following orthogonality relation: 

£ (2Z + 1)( h h 1 )( h h l ',)=SwS mm ,. (2.40) 

Using this orthogonality (2.40), we can recast expression (2.39) as 

,,2S+ 1,1,1 jsEEE (-D 



/ 7 ; 7 \ 2 

(2Zi + 1)(2Zl + 1) (2Z + l)A i_1 (Zi, S)X(Il, S - J)X 2 (l, J) 



h II I 




The explicit value of ( ^ ^ ^ ) is given by 







(2.42) 



(-1)* 



(2 5 -2i 1 )!(2 5 -2/ L )!(2 fl -20! 



(2 5 +l)! (p-Z 1 )!( 5 -Zz,)!(ff-0!' 



if Zi + Zl + Z = 2g (<? e N), otherwise zero. Finally, normalization of degenerate 
VBS states |VBS L (J,M)) is obtained as 

(VBSl(J,M)|VBS l (J,M)) (2.43) 

(2J+1)!((2S + !)!)*■ V V V 

(S + + J + l)!(5- J + l)!(J + l)!(J + l)!^ o ^ o ^ 

(2Zi + 1)(2Z L + 1)(2Z + l)A i_1 (Zi, S)A(Z L , S - J)A 2 (Z, J) 



Zi Zl Z 


Combining results of (2.21) and (2.43), we arrive at the following theorem on 
eigenvalues: 

Theorem 2.2. Eigenvalues A (J) ( J = 0, . . . , S) of the density matrix are indepen- 
dent of Cl and/or M in defining eigenvectors (see (2.15) and (2.16)). An explicit 
expression is given by the following triple sum 



A(J) 

5 + 1 
(25+1)! 



S\S\ 

j^(VBS L (J,M)\VBS L (J,M)) (2.44) 



(2j + iy.s\s\ ^ V V 

(s+j+ms-j+mj+mj+iy-tothk 



^ 2 

\L-1, 



•(2Zi + 1)(2Z L + 1)(2Z + l)X L - L (h,S)X(l L , S - J)X Z {1, J) 



Zi h I 
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Although not straightforward to verify, this expression (2.44) should be con- 
sistent with eigenvalues given through the recurrence expression (2.33). We could 
check the case when 5=1 that 

(VBS L (0,0)|VBS L (0,0)) = i(3 L + 3(-l) L ), 
(VBSi(l,M)|VBS L (l,M)) - i(3 L - (-1) L ), (2.45) 

where we have used the selection rule of the Wigner 3j-symbol. From (A 5) we find 
that (VBS|VBS) = 2 • 3 L , so that we obtain the correct eigenvalues of the density 
matrix from the above result (2.44) (see §4 for comparison). 

We shall emphasize at this point that given eigenvalues (2.44), both von Neu- 
mann entropy 

s 

S V . N = -Tr [p L In p L ] = - ( 2 J + 1) A (^) ln A (^0 ( 2 - 4 6) 
and Renyi entropy 

Sr = j^— In {Tr [pi]} = - J- In {^(2 J + l)A a (J)\ (2.47) 
can be expressed directly. 



3. Density Matrix in the Large Block Limit 

In the limit L — > oo, that is when the size of the block becomes large, we learned 
from Fan, et al. (2004), Katsura, et al. (2007a), Vidal, et al. (2003), Hadley (2008) 
that the von Neumann entropy reaches the saturated value S v .n = ln.(S+l) . 
Then the density matrix (denoted by p^ in the limit) can only take the form (see 
Nielsen & Chuang 2000 for a general proof) 

Poo = jsTT)^ I(s+1)2 ® $co ' (3,1) 

where Irs+i) 2 is the identity of dimension (S+ 1) 2 and is an infinite dimensional 
matrix with only zero entries. In this section, we give a proof of Conjecture 1 (2.30) 
in the limiting case as L — > 00. Then we shall verify the structure of the density 
matrix (3.1) explicitly. 

We first realize from (2.25) that as L — > 00, A i_1 (Z,5) — > 5^0- Therefore only 
the first term with I = is left in (2.28). So that we need only to calculate Kl(Cl): 

dtiid&L (u ia \ + (uu{ + vvl) J 

■ {u\v* L - v{u* L ) (uu* L + vv* L ) {u L al+v L blJ . (3.2) 

It is useful to know transformation properties of the integrand in (3.2) under 
SU(2). The pair of variables (u,v) defined in (2.3) and bosonic annihilation opera- 
tors (a, b) in the Schwinger representation both transform as spinors under SU(2). 
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That is to say, if we take an arbitrary element D £ SU(2) (2x2 matrix), then 
(u,v), etc. transform according to 

:H(")' 

On the other hand, (u*,v*), (—v,u), (a^,tf) and (—6, a) transform conjugately to 
(u,v). That is to say (u*,v*), etc. transform according to 

The combinations appeared in Kq(Q) (3.2) 

uia\+vib\, uui+vvl, u\v* L — ViU* L , uu* L + vv* L , uj^a\-\-vjJ^ L (3.5) 

as well as in (2.14), boundary operator in (2.11), etc. all transform covari- 
antly under SU(2), i.e. those expressions keep their form in the new (transformed) 
coordinates. 

These transformation properties (3.3), (3.4) can be used to simplify the Kq(Ci) 
integral. We first make a SU(2) transform 

D UL = ( U l !M, D U J UL ) = (l), (3.6) 



v L u L ) L \ v L J V 

under the part of the integral (3.2) over Cli. Then this part of integral becomes 

dOi (u ia \ + Vibty {uul + vv*) J {-v* 1 f~ J . (3.7) 

This can be calculated using binomial expansion and the result is 

J?-(uai+vb\y (-b\) S ~ J . (3.8) 

Then we make an inverse transform in (3.8) using = D\ L , consequently (3.2) 
is put in a form with a single integral over &l remaining: 

Kim = ^-j(ua\+vb\y (3.9) 

• J dCl L (a\v* L - b\ulj S J (uu* L + w* L ) J (u L a[ + v L bty . 



Now we make another SU(2) transform using 

"*=(-»";). M:)=(;> <3io) 

then the remaining integral over CIl in (3.9) becomes 

J dh L (a\v* L - b\u* L ) (u* L ) J (u L a\ + v L b[) * . (3.11) 
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Using again binomial expansion, this integral (3.11) yields 




(3.12) 



At last we make an inverse transform in (3.12) using D u 1 = and plug the result 
into (3.9), the final form is 

^ (A) = (sTi) 2 ^- (3 ' 13) 

This expression is consistent with Conjecture 1 (2.30), which also proves that 
{|G; J, Q}} is a set of eigenvectors of the density matrix as L — > oo. Let's denote 
the density matrix in the limit by p^. Then (3.13) leads to the result (see (2.32)) 

Poc \G;J,n) = j^-^\G;J,n). (3.14) 

We find from (3.14) that the limiting eigenvalue Aoo = rgjip is indepen- 
dent of J. Any vector of the (S + l) 2 -dimensional subspace spanned by the set 
{|G; J, ft)} is an eigenvector of p^ with the same eigenvalue rgqrrp ■ Therefore 
p^ acts on this subspace as (proportional to) the identity I( S+1 ^2. So that we 
have proved explicitly that the density matrix takes the form (3.1) in the large 
block limit. In addition, we also derive from the eigenvalues that the von Neu- 
mann entropy S v .n — — X!j=o(^ + l)Aoo In A^ coincides with the Renyi entropy 
Sr = In |Sj=o( 2J + i)^} and is cc l ual t0 tnc saturated value ln(5 + l) 2 . 

4. Density Matrix for Spin S — 1 

In the case of spin S = 1, we could prove Conjectures 1 (2.30) for finite block 
by calculating K\(£i) defined in (2.29) using similar methods as been used in §3. 
However, in this special case S = 1, we have an alternative algebraic proof. We 
shall use a different representation in which the eigenvectors of the density matrix 
form an orthogonal basis (maximally entangled states) . The formulation is base on 
Fan, et al. (2004). 

(a) Ground State of the Unique Hamiltonian 

The unique Hamiltonian is given by (1.1). In order to represent the unique 
ground state, we first introduce the following notation for convenience (Fan, et al. 
2004): 

\a) = (-l) 1+5 ='°/®a Q |0), a -0,1, 2, 3 (4.1) 

where ao = I (2-dimensional identity), cr a =i t 2,3 are Pauli matrices and |0) = ^=(| Tl 
) ~ I it)) is the singlet state (antisymmetric projection) of two spin-l/2's. The four 
states (4.1) (maximally entangled states) form an orthonormal basis of the Hilbcrt 
space of two spin- 1/2 operators. 

The spin-1 state at each site is represented by a symmetric projection of two 
spin-1/2 states given by (4.1) for a = 1,2,3. Let's take the jth site for example. 
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The two spin-l/2's are labeled by (from left to right, respectively). Then the 

spin-1 states are prepared by projecting these two spin-l/2's (4-dimensional space) 
onto a symmetric 3-dimensional subspace spanned by 



1 

-1 



1% = -^(1 Thl T>? - I i>il i>j), 



1% = -s(IT>ilT>5 + U>iU>3), 



1% = ^=(1 T)jl + \ I) A T)j)- (4.2) 

Thus the two ending spin-l/2's are labeled as site and N+ 1. The unique ground 
state in this representation is (Fan, et al. 2004; Affleck, et al. 1987, 1988) 

|G> - « =1 P,j) |0) 0l |0) l2 • • • |0)^ +1 . (4.3) 

Here Pfj projects two spin-1/2 states onto a symmetric subspace, which describes 
spin-1. Using basis (4.1), we have 

3 

Pfi = Y,\ a )fM- ( 4 - 4 ) 

a=l 

A crucial step (see Fan, et al. 2004) is that the ground state (4.3) can be expressed 
in a different form using 



-1 3 

\0)ab\0)bc = — E \ u )bb [I A ® Mc] \°)ac (4-5) 



d = 



for arbitrary labels A, B and C. Repeatedly using relation (4.5), the product of 
|0)'s in (4.3) can be rewritten as 

|0)oi|0)i2---|0)^ +1 (4.6) 

r) E hi) • • • I«jv) [h® {v aN N+1 ] |o)ojv+i- 



Then by projecting onto symmetric subspace spanned by \a — 1,2,3), the ground 
state (4.3) takes the form (Verstraete, et al. 2004a; Fannes, et al. 1992) 

1 3 

l G ) = E l«i> ■ ■ ■ l"JV> [lo®(° aN N+1 ] |0) jv+i- (4-7) 

Note that this ground state (4.7) is normalized and we have re-written the overall 
phase for it has no physical content. 



(b) Density Matrix of a Block of Bulk Spins 

Given the ground state in the form (4.7), we obtain the density matrix of a 
block of L contiguous spins starting at site k by tracing out spin degrees of freedom 
outside the block using basis (4.1): 

Pl = Tr Bili ... ifc _ lifc+L ,... iJViJV+1 |G)(G|. (4.8) 
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The result is independent of the starting site k and the total length N (see Fan, et 
al. 2004). We choose k = 1, N = L so that the density matrix reads (Fan, et al. 
2004) 

1 3 

Pl = ^l \ a i}( a i\---\ a L)(a / L\(Q\ I ® (<X ■■■o a > L )I® (v aL ■••o- ai )|0). 

a,a' = l 

(4.9) 

(c) Ground States of the Block Hamiltonian 

The degenerate Hamiltonian is given by (1.3). We choose the length of the spin 
chain to be equal to that of the block, then the block Hamiltonian Hb = Hdeg with 
N = L reads 




(4.10) 



Any linear combination of states of the following form 

|G;xi,Xz> = (®£=i^j) IXi)i|0)i 2 |0)23 ■ ■ ■ \0)— 1l \xl)l (4.11) 

is a ground state of the block Hamiltonian (4.10). In (4.11) we have made notation 
|x) = | | or J.) represents the two spin-1/2 states and Pjj is defined in (4.4). 
Let's make a particular linear combination of these \G;xi,Xl) states using (4.1) 
and write the four (a = 0,1,2,3) linearly independent ground states of the block 
Hamiltonian (4.10) as follows 

\G;a) (<8>f =1 P n ) \a) Ll \0) l2 \0h 3 - ■ -\0)— 1L . (4.12) 

Now we go through the same steps as from (4.3) to (4.7), the resultant form of the 
four ground states (a = 0, 1, 2, 3) is 

3 

|G;a}= K)---K> {aL\tr a <»(tT aL _ 1 ---a ai )\0). (4.13) 

,ot L = l 

Again we have re-written the overall phase for simplicity. These four states are 
orthogonal, and the normalization is given by 

f l(3 L + 3(-l) L ), a = 0- 
(G;a\G;a) = \ (4.14) 
{ \(3 L -(-l) L ), a = 1,2,3. 

(d) Eigenvectors of the Density Matrix 

According to theorem 2.1, the degenerate ground states (4.13) are eigenvectors 
of the density matrix (4.9). Let's apply p L to \G;a) and use orthogonality of the 
| a) states. Then we obtain 

1 3 

Pi|G;a) = ^ l"i> ■■■!«£> C ai ... aL (4.15) 
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with coefficient 

3 



C ai ... aL = KK® (o- qL i • ••o- Qi )|0) (4.16) 

•(0|J® (<7 a / • ■ -o-^)/® ( )|0). 



It can be shown by induction that 

3 



53 (/®a ai _ 1 --- f 7 ai )|0)(0|(7® ( 7 ar --a ai _ i ) = 53^|/3}(/3| (4.17) 

with 

r ^"^(-l) 1 - 1 ), = 0; 
^ = ^ (4.18) 
[ 1(3^-1 - (-l)^-i), = 1,2,3. 

Therefore the coefficient C Ql ... ai/ defined in (4.16) can be simplified as 

C ai ... aL (4.19) 

3 

53 A p (a / L \a a (g)I\f3)(P\I(g>(cr < a aL )I(g)(a aL _ 1 ---<J ai )\0). 

a> L =l,P=0 

Straightforward calculation using multiplication rules of Pauli matrices shows that 
(4.19) can be further simplified as 

$afl{(XL\I ® (ffoL-i ' ' ' °"ai)|0) (4.20) 
3 

+ (A + 2Ar)(l - <5 a , )(<W (0| - i 2 e aai/J <0|)J ® ( )|o) 

where e QQL( 3 is the totally antisymmetric tensor of three indices with £123 = 1. By 
realizing that 

3 

<W(0| -i53 e ««i-/3(' 9 l = (°\ (J aLV a ® I = (a L \<j a ®I, (4.21) 
0=1 

we have reached the final form of the coefficient C ai ... aL such that 

C ai ... aL = [3A 1 5 a ,o + {Ao + 2A 1 )(l-5 a ,Q)](a L \<T a <2)(a aL _ 1 ---a ai )\0). (4.22) 
As a result, we plug (4.22) into (4.15) and find that 

pJG;«> = 3Al ^ + (Ao+ L 2 - 4l)(1 -^ ) (4.23) 

O 

3 

53 \ai)---\a L )(a L \a a (g)(a aL _ 1 ---(T ai )\0). 

<Xl,— , a L = l 
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By comparing with (4.13), we find that (4.23) is exactly the statement that \G;a) 
(a = 0, 1, 2, 3) are eigenvectors of the density matrix p L : 

p L \G:a) = A a \G;a), a = 0,1, 2, 3 (4.24) 

with eigenvalues 

A a = 3^i^,o + (A) + 2^i)(l-^,o) = M (1 + 3 H )L) ' a = ° ; (4 25) 

3L I a = 1,2,3. 

These numbers obtained in (4.25) are exactly the eigenvalues found in Fan, et 
al. (2004), Katsura, et al. (2007a) for spin-1, and are consistent with our explicit 
expression for eigenvalues (2.44). 

We can also prove explicitly that any other eigenvectors of p L orthogonal to the 
set {|G;a)} have zero eigenvalue. Let's note that a complete basis of the Hilbert 
space Hl of the block of spins can be chosen as 

{\ ai )---\a L )}, a =1,2, 3. (4.26) 

The subspace H\ with non-zero eigenvalues is panned by {|G;a)} ; as we have 
already shown. The Hilbert space can be reduced into a direct sum 

Hl = H\ ® (4.27) 

We will show that the subspace J?$ orthogonal to H^ is a subspace of van- 
ishing eigenvalues. Mathematically, this means that for an arbitrary basis vector 
\0i) ' ' ' \Pl), we shall have 

P l {Il-Pa)\0i)---\0l)=O, (4.28) 

where I L is the identity of H L and Pa is the projector onto if a : 

Il= E K>---K)(«i|-->l|, fA^E /rTI ( 429 ) 

By taking expressions (4.9), (4.29), (4.24), and realizing that 

a=0 a=0 

we find the left hand side of (4.28) being equal to 

Pl (Il-Pa)\Pi)---\Pl) (4.31) 
1 3 

= ? ^ K}|---K) (0|[/<8> (o-ft ■••<t /3 J,/® (a aL ■••<r ai )]|0). 



ai-.-CKz,— 1 



We use multiplication rules of Pauli matrices to write the two terms within the 
commutator in (4.31) as 

I ® (<J 01 ■■■<J 0L ) = e ie(/3) 1 0(7/3, /3 = 0, 1, 2, 3; 

J® (<7 Qi •••<r Ql ) =e ie ^/®CT Q , a = 0,1, 2, 3. (4.32) 
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Here e 10 ^ and e 10 ^ are two phase factors. Then the commutator is 

[I ® (<t 01 ■ • • <7 fo ), / ® (a aL ■■■a ai )]= c^+^I ® [ CT/3 , <7 Q ] . (4.33) 



There are two possibilities: (i) a = (3 or at least one of the two is equal to zero, 
then erg and a a commutes; (ii) a ^ /J ^ 0, then [up, a a ] = 21^ , e/j ai (J 1 , but we 
still have (0\I ® ct 7 |0) = (0|7) = 0. Therefore, the factor (0|[7 ® (a^ ■ ■ ■ a 0L ), I ® 
{&a L ■ ■ ■ °ai)]|0) m (4-31) is identically zero. So that we have proved (4.28). There- 
fore H<j> is a subspace with only zero eigenvalues. 



5. A Different Proof of the Theorem on Eigenvectors 

It was shown in §2 d that the density matrix takes a diagonal form in the basis of 
zero-energy ground states of the block Hamiltonian (2.13). In this section, we show 
the same result by taking a different approach. This alternative proof of theorem 
2.1 does not involve coherent state representation. 

Let's start with the ground state of the unique Hamiltonian (1.4) with N = L: 



|VBS) = J] (a]b] +1 -6].a]. +1 ) S |vac). (5.1) 
3=0 



In order to calculate the density matrix p L = TY ,l+i [p], where p is defined in 
(2.2), we introduce a useful identity: 



(_l)S-J+M 

o,l + i(J,M\(\8) ,i®\8)l,l+i) = + \ J >- M Kl, (5.2) 



where \J,M)o,l+i is identical to the spin state defined in (B6) except for site 
indices. \s)ij in (5.2) is the normalized singlet state with S valence bonds defined 
as 



= A=£= ]T (-inSA-m^lSAm^. (5.3) 

V 5 + 1 m =-S/2 



Identity (5.2) is derived using properties of the singlet state (5.3) and Clebsch- 
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Gordan coefficients as follows: 

o,L+i(J,M\s) aA \s) L ^ L+1 

m a +m L + 1 =M 

]T (J,M|S/2,m ;S/2,m L+1 ) (S/2 

m ,m L + 1 

f-n- s/2 

■ -^L J2 (-iri^A-mOo^Am!)! 

mi = — S/2 

r-n- s/2 

• A=£= J] (-1)^15/2, -m L ) L |5/2,m L ) L+1 

m +m L + i = M 

= E (-l) mo+mi - +1 (^M|5/2,mo;5/2,m L+1 ) 

• |5/2,-m )i|5/2,-m i+ i) L . (5.4) 
Here the Clebsch-Gordan coefficient is defined by 

(J, M\S/2, m ; S/2, m L +i) = ij(J, M\ {\S/2, m ), ® \S/2, m L+1 ) 3 ) . (5.5) 
Then using the symmetry property of Clebsch-Gordan coefficients 

(J,M\S/2,m -S/2,m L+1 ) = (-1) S - J ( J, -M \S/2, -m ; S/2, -m L+1 ), 

(5.6) 

and the completeness of the basis {|S/2,m )o ® |S/2, m L+1 ) i+1 }, we obtain the 
identity (5.2). 

With the help of identity (5.2), we calculate the partial inner product of the 
VBS state with the state | J, M)o,l+i, which is involved in taking trace of boundary 
spins. The VBS state |VBS) is decomposed into the bulk part and edge parts, then 
making use of (5.2), we have 



o,l+i(J,M|VBS) 

L g 

= 0>L+1 (J,M\]l(a]b] +1 -b]a] +1 ) |vac) 



3=0 
L-1 



3=1 

L-1 



= J] (a]b] +1 -b]a] +1 ) o,i+i<^M|*)o,i|«)i,i + i|vac) 2 ...i-i 



= ( 5! ) 2 LI ( a >U b Wj + i) (-1) S - J+M \J, -M) 1 , i |vac) 2 ... L _ 1 

3=1 

= (-1) S - J+M (S!) 2 |VBS L (J,-M)). (5.7) 

We see that the (S + l) 2 degenerate VBS states |VBS L (J,M)) defined in (2.16) 
appear in the partial inner product (5.7). As discussed in §2 c , they form a complete 
set of zero-energy ground states of the block Hamiltonian (2.13). These states are 
nothing but the edge states of the subsystem (block). 
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Now, it is straightforward to evaluate density matrix as 



Tr .L+i[p] = Yl 




o,l+i{J, M|VBS) (VBSL7, M)q, l+ i 
(VBS|VBS) 



^|VBS L (J,-M))(VBS L (J,-M)|. (5.8) 



J,M 



This expression is identical to (2.21) as we change dummy index from M to —M. 
Therefore, in this approach again we arrive at theorem 2.1 that the density matrix 
is proportional to a projector onto a subspace spanned by the (S + l) 2 ground states 
of the block Hamiltonian (2.13). Normalization (VBS|VBS) is given in Appendix A. 
States |VBSl(J, M)) are shown to be mutually orthogonal in Appendix C. 



We have studied the density matrix p L of a block of L contiguous bulk spins in 
the AKLT model. The unique Hamiltonian for generic spin-S 1 is given by (1.4), 
which has a unique ground state described by the VBS state (2.1) in the Schwinger 
representation. The density matrix p L (2.10) of the block is obtained by taking 
trace (2.6) of all spin degrees of freedom outside the block. The structure of the 
density matrix has been investigated both for finite and infinite blocks. 

For generic spin-S 1 and finite block, two mathematically rigorous results have 
been established as theorem 2.1 and 2.2. In theorem 2.1 we constructed eigenvec- 
tors of the density matrix with non-zero eigenvalues. These eigenvectors |G; J, O) 
defined in (2.15), or |VBSz,(J, M)) defined in (2.16) equivalently, are proved to be 
the (S + l) 2 zero-energy ground states of the block Hamiltonian (2.13). The cor- 
responding eigenvalues are obtained in two different forms. Using nonorthogonal 
basis {|G; J, fi)}, the eigenvalues are given through Conjecture 1 (2.30) and the re- 
currence relation (2.31); while using orthogonal basis {|VBSl( J, M))}, in theorem 
2.2 an explicit expression (2.44) for eigenvalues in terms of Wigncr 3j-symbols is 
derived. Non-zero eigenvalues A (J) with J = 0, 1, . . . , S ((2.33) and (2.44)) depend 
only on J and are independent of f2 and/or M in defining eigenvectors. The density 
matrix (2.21) is a projector onto the subspace of dimension (S + l) 2 spanned by 
the set of eigenvectors {|G; J,Q}} and/or {\VBS L (J, M))}. 

In the large block limit L — > oo, Conjecture 1 (2.30) is proved and all non-zero 
eigenvalues Aoo become the same (3.14). The infinite dimensional density matrix p M 
(3.1) is a projector onto a (S + l) 2 -dimensional subspace in which it is proportional 
to the identity. The von Neumann entropy S v .n coincides with the Renyi entropy 
Sn and is equal to the saturated value ln(5 + l) 2 . In the limit the Renyi entropy 
is a independent, which behaves quite differently from the XY model where the 
Renyi entropy has an essential singularity as a function of a (see Its, et al. 2005; 
Franchini, et al. 2007, 2008). 

We have also investigated the structure of the density matrix in a special case 
when spin S — 1. Both theorem 2.1 and 2.2 are proved using a different represen- 
tation (4.13) (maximally entangled states) where all four eigenvectors \G,a) are 
orthogonal. We have also shown (4.28) explicitly that any vector orthogonal to the 
subspace spanned by the set {|G, a)} has zero eigenvalue. 



6. Conclusion 
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Based on the main results obtained in this paper, we end our conclusion by 
making the following conjecture: 
Conjecture 2 

The structure of the density matrix as a projector onto a subspace is gener- 
alizable to inhomogeneous AKLT spin chains (spin values at different lattice sites 
could be different) and lattices of higher dimcnsionsf . In the large block limit, the 
density matrix should behave as the identity operator within the subspace. i.e. Eq. 
(1.7) is valid for arbitrary large lattices. 

The authors would like to thank Professor Heng Fan, Professor Anatol N. Kirillov and 
Professor Sergey Bravyi for valuable discussions and suggestions. The work is supported 
by NSF Grant DMS-0503712 and the Japan Society for the Promotion of Science. 



Appendix A. Normalization of the VBS State 

The VBS state |VBS) (also known to be the ground state of the unique Hamiltonian 
(1.4)) defined in (2.1) is not normalized. Using the coherent state formalism (2.4) 
and the completeness relation (2.5), we express the norm square as 



(VBS|VBS) 



(5 + 1)! 



47T 



where we have used 



2 r 



(2 5+1)! 

47T 



N 



r N + l \ N 

/ n<^ n 



; (i-n r si j+1 ) 



(0\a s+l b s - l \Cl) = ^/{2Sy.u s+l v s - 1 . 



(Al) 



(A 2) 



(A3) 



Now we expand i(l — Clj ■ Clj+i) in terms of spherical harmonics as in (2.24), 

then integrate from tl to fijv+i- We notice by using the orthogonality of spherical 
harmonics that each integral contributes a factor of except the last one. For 
example, 



/ 



d^ 

47T 



j(i-fio-ni) 



s + 

47T 

5 + 1 



T J>(Z,5) ]T v^Y^i) / dn Y lm (Q )Y* Q (n Q ) 

1=0 m=-l J 



AttY*^ 



Air 
5+1' 



(A 4) 



The last integral over Oat+i contributes simply a factor of 47r. Consequently, the 
norm square (A 1) is equal to 



(VBS | VBS) 



(25+1)! 



5!(5 + 1)!. 



(A 5) 



f It may even be gencralizablc to certain classes of arbitrary graphs if the limiting sub-graph 
can be defined properly. 
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Appendix B. Rank of the Set {|G; J, Cl)} with Fixed J Value 

For notational convenience, we define 



J+M..J-M ( r A\S+mfiA\S-m 



*L= 7=== ., (Bl) 



y/(J + M)\(J-M)l' bm y/(S + m)\{S-m)\ 



These two variables transform conjugately with respect to one another under SU(2). 
Xjm has the following orthogonality relation 



/ 



47T 

dflX} M Xj M ' = (2J + iy &MM' ■ (B 2) 



ipg m is a spin state creation operator such that 

4jvac) = |S,m>. (B3) 
The operator A\ defined in (2.14) can be expanded as (see Hamermesh 1989) 



4 = ^^'^-O^: (B4) 



(5 + J + l)!(5- J)\J\J\ 
2J + 1 

^ J! (5/2,mi;S/2,m2|J,M)^ Arni ®4 /2imj ., 



M=— J mi,mi 

where (5/2, mi; 5/2, m 2 | J, M) are the Clebsch-Gordan coefficients. Note that 

mi anc " mi are defined in the Hilbert spaces of spins at site 1 and site 
L, respectively. We realize that the particular form of the sum over mi and rriL in 
(B 4) can be identified as a single spin state creation operator 

^jm= E {S/2,m 1 ;S/2,m 2 \J,M) ^ s/2mi ®i>s/2,m L - ( B5 ) 

mi ( mt 

This operator acts on the direct product of two Hilbert spaces of spins at site 
1 and site L. It has the property that 

* t 7M |vac)i®|vac) L = |J,M) liL . (B6) 

Now we can derive the completeness relation of the set {|G; J, tl)} using (B 2), (B 4) 
and (B5): 

J dCl\G; J,Cl)(G; J,Cl\ (B7) 
4tt (5 + J+l)!(5-J)!J!J! ' t 

- WTTy. 2TT1 1,* jm \vbs l )(vbs l \*j M . 

The set of states {4 , j M |VBS.l), M = — J,..., J} are linearly independent. So 
that the rank of {|G; J, 0)} with fixed J value is 2 J + 1. With the introduction of 
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degenerate VBS states |VBSl(J, M)) in (2.16), |G; J, Q) can be written as a linear 
superposition: 

M=-J 

More details can be found in Hamermesh (1989). 

Appendix C. Orthogonality of Degenerate VBS States 

The set of degenerate VBS states {|VBS L ( J, M)>, J = 0, S, M = -J,..., J} 
introduced in (2.16) are mutually orthogonal. To show this, it is convenient to 
introduce the total spin operators of the subsystem: 

L L L 

s+ t = X>]^ Stot = E 6 k-< = ^>k- - h h)i 2 - ( C1 ) 

i=i i =1 o =1 

First we show that the set of operators {S^~ ot , S^ ot ,S^ ot } commute with the product 
of valence bonds, i.e. 

n>5&5+i - 6 5°5+i) s ] = °> t^tot, nV^+i - b hUf] = °- ( c 2 ) 

i=i i=i 

These commutation relations (C 2) can be shown in similar ways. Take the commu- 
tator with S^ ot first. We re-write the commutator as 

= X>M - b\4) s ■ ■ ■ [S+ , (a]b] +1 - b]a] +1 ) s ] ■ ■ ■ (a^bl b L _A) S 
i=i 

= X>t4 - &t4) S • • • ^ + (aj6t +1 - b]a] +1 ) s ] . . . 



■■(4-i^-6Li4) f 



Then using commutators [dj, a|] = % and [6j, bj] = <5y, we find that 



(C3) 



Wjbj + a] +1 b j+1 , {a)b) +1 b)a) +l ) s ] 

at(-S)aJ +1 (aj6t +1 - &J a t +i) s-i + a t +i5a t (a t 6 t +i _ fcj t +i) 5-i 

0. (C4) 
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Therefore [S? ot ,U*f=i( a ] b ]+i - b ] a ]+i) S ] = °- In ( C4 ) wc have used i b 3> Aj b ]+i ~ 
b l a l+i) S ] = — ^ a ]+iAj b l+i ~ kjaj+i) 5-1 - ^ n a P ara Uel way, we find that the com- 
mutator with 5 t ~ t also vanishes. Next we consider the commutator with S* ot : 

{SLul[( a l b U- b H +1 ) s } 

3 = 1 

= X>i4 - b \4) S ■■■[S] + S] +1 , (a]b] +1 b]a] +l ) s ] ■ ■ ■ 

■■■^L-A-A-Af- 

(C5) 

In the right hand side of (C 5) , the commutator involved also vanishes because 

[s> + s> +l M) b ] +1 - b A + A\ 

= -j[ a ] a 3 - b ] b 3 + a ) + l a 3 + l - b j+l b j+U AA3 + I ~ b A+A] 

= a ][ aj , AA+i b A+A} - AA+i b AA s ] 

+ a] +1 [a 3+u (atftt +1 - b]a] +l ) s ] b] +1 [b j+1 , {a)b] +1 b)a\ +1 ) s ] 

= (C6) 

Substituting (C6) into (C5), wc obtain [S^ ot , Uj=i( a ] b ]+i ~ = °- Now 

we shall show that the state |VBSl(J, M)) is an eigenstate of S* ot and the square 
of the total spin S 2 tot = ^(S^ ot S^ ot + S^ t S^ t ) + (5 t z ot ) 2 with eigenvalues M and 
J(J + 1), respectively. Using the commutation relations (C2), we can show that 



L-l 



S± |VBSi(J, M)) = [] (a]b] +1 - b]a] +1 ) s (Sf + S±)\J, M} 1 , i |vac} 2 
3=1 

L-l 

SZJVBS L (J,M)) = H(a]b] +1 - bla] +1 ) s (S z 1+ Sl)\J,M} hL \v a c) 2 _ 



L-l 



L-l- 



3 = 1 



(C7) 

Then from the definition of the state |VBSl(J, M)) and the following relations: 



(S+ + S+)\J,M) hL = y/(J T M)(J ± M + 1)| J, M ± 1), 

(S z + S z L )\J,M) hL ) = M\J,M) hL , (C8) 

we obtain 

S£t\VBS L (J,M)) = y/(J tM){J±M + l)|VBSi(J, M ± 1)), 

S t z ot |VBS L (J,M)) = M\VBS L (J,M)) (C9) 

and hence S^ ot \VBS L (J, M)) = J (J + 1)|VBS L ( J, M)>. It is now proved that 
|VBSl(J, M)) is an eigenstate of S* ot and Sj ot with eigenvalues M and J(J + 1), 
respectively. Therefore the states with different eigenvalues (J, M) are orthogonal 
to each other. 
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Appendix D. Density Matrix and Correlation Functions 

The relation between the density matrix and correlation functions was sudied in 
Jin & Korepin (2004), Katsura, et al (2007a), Arovas, et al (1988). It was shown 
in Section 2 of Jin & Korepin (2004) that the density matrix contains information 
of all correlation functions in the ground state. The original proof was for spin 
5 = 1/2. In this appendix we shall generalize this result to generic spin-5 which is 
applicable to our AKLT model. 

The Hilbert space associated with a spin-5 is (25 + l)-dimensional. Therefore 
we could choose a basis of (25 + l) 2 linearly independent matrices such that an 
arbitrary operator defined in the Hilbert space can be written as a superposition 
over the basis. Let's denote the basis by {A ab ; a, b = 1, . . . , 25 + 1}, in which each 
matrix A ab is labeled by a pair of indices a and b with totally (25 + l) 2 possible 
combinations. The matrix element is defined as 

(A ab ) M = S ak S bl , k,l = l,...,2S+l. (Dl) 

In addition to {A a h}, we introduce an equivalent "conjugate" basis {A a &} such that 

{A ab ) k i = 6 a i5 bk , a,b,k,l = 1,...,25 + 1. (D2) 

These matrices (D 1) and (D2) are actually matrix representation of operators 
{|5, to) (5, m'\; to, to' = —5, . . . , 5}. They arc normalized such that 

Tr(A ab A cd ) = ^(A ab ) k i(A cd )i k = ^S a iS bk 5 c iS dk = S ac 5 bd . (D3) 
k,l k,l 

Here TV takes trace at one and the same site. Because of the completeness of {A a &} 
at each site, the density matrix of the block can be written as (see (4.8)) 

Pwock = 7>outside|G}(G| = (®?e {Wock} A aj b] ) coeE{a 3 b }, (D4) 

{a jb] } 

where |G) denotes the unique ground state, Tr ou tside takes traces of sites outside 
the block and coefflajOj} denotes the coefficient. Using the normalization property 
(D 3), the coefficient coeff{aj&j} with label j taking values within the block can be 
expressed as 

coeS{ aj bj} = Y II Tr{A a]b] A C]d] )cocR{c 3 d 3 } 

{cjdj} jGblock 
= TVblc-ck [ (deblock Aj^) Pblock] 
= Tr a ll [(® J - eb lockAy&,) |G)(G|] 

= (G| (<8>jgbiock4.,6,) |G). (D5) 

Here Tr b i ock takes traces of sites within the block and Tr a n takes traces of all lattice 
sites. Combing (D4) with (D5), we have the final form 

Pwock= (® je {block} (G| (deblock Aajbj) |G>. (D6) 

{a j bj} 

This is the expression of the density matrix with entries related to multi-point cor- 
relation functions (G| (<Sij ebiockA ajbj ) |G) in the ground state. All possible combi- 
nations {cijbj} are involved in the summation. Therefore, we have prove for generic 
spin-5 that the density matrix contains information of all correlation functions. 
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